Interaction of a point charge with the surface of a uniaxial dielectric 
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We analyze the force on a point charge moving at relativistic speeds parallel to the surface of a 
uniaxial dielectric. Two cases are examined: a lossless dielectric with no dispersion and a dielectric 
with a plasma type response. The treatment focuses on the peculiarities of the strength and direction 
of the interaction force as compared to the isotropic case. We show that a plasma type dielectric 
can, under specific conditions, repel the point charge. 



I. INTRODUCTION 

Despite the long and rich history of theoretical studies 
on the interaction between fast charges and solid sur- 
faces (see, e.g. refs. [TUl2]) unexpected results can be 
and indeed are derived. A recent example is the dis- 
covery by one of the authors of the present letter [T5] . 
that the interaction between a relativistic charge packet 
and a metal or dielectric surface can become repulsive by 
simply tuning the packet geometry; a result that seems 
to go against common notions established in electrody- 
namics and should be of importance in the framework of 
accelerator physics and electron spectroscopy. 

In this letter we switch gears and focus on the inter- 
action between a point charge and a uniaxial dielectric 
within the context of ionic and molecular interactions 
with macroscopic surfaces. We assume a description of 
the surface that approximates the non-isotropic nature 
of crystalline surfaces and surfaces decorated with ad- 
sorbed non-isotropic inclusions. We present a deriva- 
tion of the force on the charged particle starting from 
Maxwell's equations in Fourier space and then evaluate 
the force numerically in real space. We show that the 
longitudinal component of the force (parallel to the di- 
electric surface) is in general not parallel to the particle 
velocity anymore and that its direction depends on the 
particle speed (energy). We demonstrate two peculiari- 
ties of the plasma type response: 1) the direction of the 
longitudinal force depends on the distance of the parti- 
cle from the surface and 2) under specific conditions the 
particle can be repelled by the surface. 



II. EVALUATION OF THE 
ELECTROMAGNETIC FORCE 

The geometry of the problem is illustrated in fig. [T] 
A point charge moves in vacuum with a velocity v at a 
distance zq parallel to the surface of a uniaxial dielectric. 
The dielectric surface lies in the xy plane with the op- 
tical axis oriented along the x direction. The velocity is 
V = {vx,Vy,Qi) = (ucos0,usin0, 0), where 9 is the angle 



between v and the optical axis. 

The electromagnetic field due to the moving point 
charge is calculated from Maxwell's equations: 



VD = p 



V X H = J 






(1) 
(2) 



The problem is tackled by replacing the fields with 
the standard scalar and vector potentials defined as B = 
V X A and E == -V$ - dA/dt. Then the solution to 
eqs. (IT]) and (pi) is sought separately in the vacuum (z > 
0) and dielectric {z < 0) half spaces by introducing three- 
dimensional Fourier transforms of all quantities: 



G(r,z,i) 



cPkdu} 

(27r)3 



g(k,z,a;)e'(''- 



(3) 



where k = (/c^,, ky) and r = {x,y) are the wave and posi- 
tion vectors parallel to the dielectric interface. 

The absence of bound charges and currents in vac- 
uum allows us to decouple the above equations using the 
Lorenz gauge: 



a2$i 



dz 
dz"^ 



2 -Ql$l = 

-g?Ai = 



J 

c^eo' 



(4) 
(5) 



where Qi = -^/fc^ „ uS^ jc? , k = w/c^ -I- fc^ is the magni- 
tude of k and the index 1 refers to the vacuum half space. 
The Fourier decomposition of the charge density is: 



p(k, z,uj) — 2i:q5{uj — v • 'k)6{z — zq). 



(6) 



while the current density is simply 3 — pw. 

The solution to eqs. (Id]) and (Is]) is written as a sum of 
the "incident" field due to the point charge (the same as 
in the absence of the dielectric interface) : 



$i = -^5{uj - V • k)- 
eo 

V 

A,; = ^$,;, 



-Qi\z-zo\ 

~Qi 



(7) 
(8) 
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FIG. 1. The point charge moves in vacuum at a distance zo 
parallel to the surface of a uniaxial dielectric defined by the xy 
plane. The velocity makes an angle 9 with the optical axis {x- 
direction). The dielectric tensor is diagonal with e^x ~ ei('^) 



and iyy = ezz — €2(10). 



^(6iA2)fc2+fc2_,^(^2/,2) 



k 



— COS^ ( 
£2 



- sin — ei/32 cos^ 



where in- 



dex 2 refers to the dielectric half space. The above 
equations are consistent with results obtained by other 
authors (see, e.g. [HI US]), except that we prefer to work 
with potentials rather than fields. 



The general solution to eq. ( 10 ) is written as a sum of 



ordinary (o) and extraordinary (e) waves 



where 



A, = a„e'32 



ap.e 



Q2,Z 



(11) 



and the "scattered" field due to the dielectric: 



a.e 



'Qiz 



(9) 



where for the latter we impose the gauge <i>s = and 
therefore V • A^ = 0. 

In these equations Qi is a real number. This fol- 
lows from the fact that all the fields carry the prefactor 
(5(a; — V • k), which insures that w == v • k. Writing k — 
(A:cos0, ksiiKJ)) we obtain Qi = k^\ — p^ cos^ (Q — (/)), a 
real quantity for any j3 = w/c, and </>. The point charge 
moving at constant speed parallel to the dielectric sur- 
face cannot emit radiation in the vacuum half space - the 
waves in eqs. |7]) - ([9]) are evanescent and the sign in front 
of Qi in eq. (pTis negative to insure that the waves decay 
to zero as z — )• 00 

The following linear constitutive relations are assumed 
for the dielectric: D = cqC • E, where e is diagonal with 
Sxx = ei(w) and Cyy = e^z = £2(0;) and H = B//io (the 
material is non-magnetic) . From eq. (IT]) we can again set 
<I>2 = and therefore V • (e • A2) = 0. The equation for 
the vector potential in the dielectric obtained by taking 
the curl of eq. (pi) then becomes: 



d^A2z 



QleA2 



A2y - 



i2o 



A2. - 



/2o^2^ 



1 I n,x'^y^2x 



-1 i-ikx) 



dA 



2x 



dz 



(10) 



where Q20 
k^Jl — £2/52 cos 



Vfc2-e2(w7c2) 
and Q2e 



So — (0, Cloy^ 0,oz) 



Kxf^y 



~ikxQ2 



1,2 



£2 



^2 



£2 



(12) 



Since both ei(w) and e2('^) are in general complex 
quantities, Q20 and Q2e are also complex. There are two 
solutions for Q20 and Q2e but the physical ones corre- 
spond to those with positive real parts (only these decay 
exponentially in the dielectric and satisfy the radiation 
condition) . 

To find the coefficients contained in the vectors a^, &„ 
and ag we impose the usual boundary conditions for the 
fields at the interface. The procedure, although straight- 
forward, is tedious and will not be reproduced in detail. 
Using the obtained coefficients the Lorentz force compo- 
nents are: 



fx = iqe ^^"^ [ioasx + Vyik^a^y - kyUsx)] 



fy- 


iqe 
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;he Fourier transform 
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(13) 



■" 2c2eo P' 

where R and P are defined as: 



(14) 



R = -e2(v • \<if\vlQ2o{Q2e - Ql){Q2o + Ql) + t'y(32o(Q2o - Ql){Q2e + Ql) + 2(e2 - l)fc^fcy W^Wj/Ql 



+ (£2 - l)fc; <(Q2e - Ql) + v'{Q2e + Ql) \+ c'kiQ2o{Q2o + Ol)(Ql - e202e) 



+ c^( 2(62 - l)klkyV^VyQi + e2klvl{Q2e + Qi) Q2o{Q2o - Ql) + (£2 - l)kl 



K<QiQ2c 



vl{Q2o + Ql)+vliQ2o-Ql) 



£2*320 



vl{2Q2e - Ql){Q2o + Ql) + vlQ2eiQ2o ^ Ql 



+ (£2 - l)kl 2 {-vl + vl) Ql + t2vl{Q2e + Ql) ^ + Ik^kyV^Vyi (£2 - \)kl (£2 - l)Ql 



e2Q2e +Q20 e2Ql(Q2e - Ql) +Q2o(e2Q2e - Ql) 



(15) 



P = C^klQ2o{Ql + Q2o)(Ql + e2Q2e) ^ e2(Ql + Q2e)(v • k)^ [Q2o(Qi + Q20) + ^^(£2 - 1)] 



(16) 



The force in real space is obtained by integration of 
the above expressions over k = (fccosc/), /csini/i), setting 
r — vt. For a lossless and dispersionless material the 
transform over k is performed analytically, giving an in- 
verse second power dependence of the force on the dis- 
tance zq from the surface, while the transform over has 
to be performed numerically. When dispersion and/or 
losses are included the integration can be performed only 
numerically. 

The longitudinal force fp = (/j,, fy) can be interpreted 
in Fourier space in a convenient way. The point particle 
excites electromagnetic waves in the semiinfinite dielec- 
tric and each of these waves carries a momentum propor- 
tional to k. This momentum has to be balanced by the 
particle which results in a force parallel to the interface. 
The magnitude of the momentum is determined by the 
boundary conditions and material properties. The longi- 
tudinal force Fp — {F^, Fy) in real space is then obtained 
by integration over the momenta of all the excited waves. 

In the following we treat two examples of dielectric 
response: a lossless dielectric with no dispersion and a 
plasma type response. Mathematically both cases can be 
analyzed using the Drude model for the dielectric tensor: 



£l(w) = £ool 
£2(w) = £oo2 



w: 



pi 



uj{uj + 171) 

^p2 



(17) 



where £001,002 are dielectric constants at high frequencies 
(w — >■ 00), Wpi,p2 are the plasma frequencies and 71,2 are 
the damping coefficients. For a lossless dielectric Wpi, 
a;p2 and 71, 72 are made infinitely small, i.e. £1,2(1^) = 
£ooi,oo2, while for a plasma type dielectric we take £ooi = 
£002 — 1- For reasons of consistency the frequency has to 
satisfy: 



w > 71, 2- 



(18) 



It can be shown (see, e.g. [16]) that the range of 
frequencies a point charge moving above a solid excites 
near its surface is proportional to 7w/zo, where 7 = (1 — 
v^ l(?)~^l'^ is the relativistic factor. The condition of 
eq. ( 18 1 therefore reads: 



w„ 



711 



>7l,2- 



(19) 



In addition to the above, for non-zero losses the plasma 
model is not applicable near cjpi,p2, where the imaginary 
part of the dielectric function dominates the response. 

For a lossless dielectric with no dispersion, i.e. £1(0;) = 
£ooi and £2(0;) = £oo2 are constants, Fp will be non-zero 
only if the Cerenkov condition, either for the ordinary 
or extraordinary waves (or both), is satisfied. This oc- 
curs when either Q20 or Q2e becomes imaginary, i.e. the 
waves become propagating. If both Q20 and Q2e are real, 
the waves excited in the solid are evanescent (they decay 
exponentially in the solid) and these do not contribute to 
Fp . The reason is that Fp in general decreases the energy 
of the particle and this can occur only if the charge emits 
radiation into the dielectric. 

It is straightforward to show that Q20 becomes imagi- 
nary if the particle moves faster than the critical speed: 



&o = Vo/C - 



(20) 



Above Po propagating Cerenkov waves are emitted into 
a circular cone defined by: 



1 



cosa„ 



^y^ 



(21) 



Here Uo is the angle between the optical axis and the 
wave vector K = {kx,ky,BgTi[LS\kz), and k^ — iQ2o- The 
function sgn(a;) insures that the energy fiow is directed 



into the dielectric (propagating waves are actually emit- 
ted only into one half of the Cerenkov cone) . The symme- 
try axis of the ordinary Cerenkov cone is always parallel 
to the particle velocity. 

For extraordinary waves the analysis is a bit more cum- 
bersome [IT]- The critical speed for Cerenkov emission 
is: 



1 



TeT 



sin 9 + e2 cos^ 9 



(22) 



From eq. (22) it follows that for given ei, 62, and /3 there 



exists a critical angle 9c above which the Cerenkov con- 
dition for extraordinary waves is satisfied: 



cos9r 




(23) 



The extraordinary Cerenkov cone has an elliptical cross 
section [17^ , which follows from the definition of Q2e- The 
symmetry axis of the cone lies in the xy plane and makes 
an angle x with the x axis f l8j : 



tanx 



ei - €2 + v/(ei - £2)^ + (/3^eie2)^ + W^eM^2 - £1) cos (29) 
I3^eie2sm{29) 



cot (20), 



(24) 



which for (3 — (3^ and 9 = 9c reduces to: 



ei 



tan X = — tan ( 

£2 



(25) 



For a plasma type dielectric Fp is in general always 
nonzero, but is appreciable only in a certain frequency 
range due to condition (18). This can be seen for the 
limiting case when 71 = 72 = and therefore ei_2(w) = 
1 — u}pi2/^^- If follows that Q20 is real for all w and 
there are no ordinary propagating waves. On the other 
hand Q2e is imaginary if the frequency w lies between the 
plasma frequencies [H]: e.g., tUpi < u) < uip2 for LUp2 > 
LOpi. We will refer to this frequency range the Cerenkov 
band. For zero losses, only waves in the Cerenkov band 
contribute to the longitudinal force, Fp. 
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FIG. 2. Different scenarios for the direction of the longitudi- 
nal component of the Lorentz force acting on the electron. 



III. RESULTS 

In this section we plot the magnitude of Fp, the angle 
between Fp and the electron velocity, and the transverse 
force Fz- For an isotropic material, where £1(0;) = €2(0;), 
Fp is in the direction opposite to v, while for a uniax- 
ial material this is no longer the case. We plot '0 = 
arctan(Fj,/Fj;) as a function of 9 to show the effect of 
the anisotropy (the angle between v and Fp is equal to 
tA + tt). Figure [2^ sketches two possible scenarios that may 
occur. 

In fig. [3] we show the results for a dispersionless mate- 
rial with no losses. We chose ei = 8 and t2 = 4. In the 
left column we plot -0 as a function of 9 for different values 
of /3 = v/c, while in the right column we plot the mag- 
nitude of the longitudinal force and the transverse force, 
both normalized with respect to q^ /{IGireoZg) (static im- 
age charge force for a charge above a metal surface). 

The Cerenkov condition for ordinary waves, /3o > 0.5, 
is only satisfied for bottom panels in fig. [31 For extraor- 
dinary waves eq. (23) gives 9c ~ 29°. From fig. [3J Fp is 



zero below this value and increases with 9. The direction 
of Fp strongly departs from the isotropic case and for 
9 — 9c coincides with the direction of the Cerenkov cone 
given by eq. (25). This follows directly from eq. (14 1. 
For 9c the waves are emitted only into one direction (the 
Cerenkov cone becomes a line); therefore kx and ky are 
proportional and the ratio Fy/F^ is the same as ky/kx- 
The longitudinal force is thus parallel to the symmetry 
axis of the Cerenkov cone. For 9 > 9c the direction of Fp 
departs from that of the cone. 

The transverse force i^^ is a result of a complex in- 
terplay between evanescent and Cerenkov interactions 
(see [inillS] for an explanation of the isotropic case); nev- 
ertheless, Fz only slightly varies with 9 (within 10% in the 
interval 9 G [0,90°]). As in the isotropic case ,^13], F^ is 
always attractive (negative) for a point charge; it cannot 
be made repulsive simply by increasing /3 or changing 
the dielectric constant. However, it can become repul- 
sive by replacing the point charge with a transverse line 
of charge. 



0.015 




FIG. 3. Results for a dispersionless and lossless dielectric with ei = 8 and £2 ~ 4. Left column: t/) — arctem {Fy / Fx) as a 
function of 9 for /3 — 0.45 (top), /3 = 0.5 (middle) and /3 = 0.7 (bottom). Solid black line represents the results, solid gray line 
is the isotropic case {tp = 9, valid for /3 above the Cerenkov condition) and dashed gray line is the direction x of the symmetry 
axis of the extraordinary Cerenkov cone. Right column: magnitude of the longitudinal force (solid line) and transverse force 
(dashed line) as a function of 9 ioi P = 0.45 (top), 13 = 0.5 (middle) and 13 = 0.7 (bottom). 




FIG. 4. Results for a plasma type dielectric (71.2 = 10 
components (right column). Top panels: ajp2 ~ SoJmax, < 



'^piJ 



xp — arctan (Fj^/Fj;) as a function of 9 (left column) and force 



jjpi = 10^ cjp2. Bottom panels: LUp2 = O.Scumi 
black line in the left column represents the results, while the solid gray line is the isotropic case (i/> — 1 
column is the magnitude of the longitudinal force, while dotted line represents the transverse force. 



Ljpi — 10 Ljp2. Solid 
Solid line in the right 




FIG. 5. V — axctan (Fy/Fx) (top panel), magnitude of the 
longitudinal force (middle panel), and transverse force (bot- 
tom panel) as a function of 8 for a plasma type response with 
LOpi = 0.25(jJmax, ijJp2 = 10~^a;pi. Solid black line in the 
top panel represents the calculated values, dashed line is the 
isotropic case. In the bottom panel the transverse force (thick 
gray line) is split into the Cerenkov (dashed) and evanescent 
(dotted) contributions. 



For /3 = 0.5 the Cerenkov condition for extraordinary 
waves is satisfied for all 9. For low angles "0 is above 
the value for the isotropic case: -0 > ^- At some critical 
angle we enter the regime ip < 9. Increasing leads 
to transition back to the regime ip > 0. This pecuhar 
behavior cannot be qualitatively explained by considering 
the direction of the Cerenkov cone; the direction of the 
force has to be determined by integration over all the 
momenta of the waves excited in the solid. 

For higher j3 the Cerenkov condition for ordinary waves 
is fulfilled. These waves contribute to a force in the di- 
rection opposite to the particle velocity. The direction 
of Fp therefore approaches that of the isotropic case, as 
demonstrated in the bottom panel of fig. [3j 

We note that for the two extreme points 9 = Q° and 9 — 
90° the longitudinal force always points in the direction 
opposite to V, which follows from the symmetry of the 
situation. 

As /3 increases the particle excites more Cerenkov 
waves and therefore the magnitude of Fp increases, as 



indicated in the right column of fig. [3J For /3 = 0.45 the 
maximum value of _Fp is ~ 2 orders of magnitude smaller 
than i^j,; however, the components become comparable 
at ^ = 0.7. 

For the reversed situation, ei = 4 and €2 = 8, the 
anisotropy is much less pronounced since with increasing 
/3 ordinary Cerenkov waves are excited first and these 
contribute to a force in the direction opposite to v. The 
effect (not shown here) is similar to the situation in the 
bottom panel of fig. [3) left column, except that ip < 9. 

Figures [4] and [5] show the results for a plasma type 
dielectric: 71,2 = 10~^ min(a;pi,u;p2, Wmax)- We had 
to work with finite losses to insure numerical conver- 
gence. Although the conditions for a plasma type dielec- 
tric are violated at oj = and near a;pi^p2, insuring that 
i-^max ^ 7i,2 and putting the plasma frequencies "far" 
apart (e.g., ujp2 — lOwpi), the contribution to the inte- 
gral around these points is small (the numerical results 
are nevertheless exact because we are using the Drude 
model for integration). 

Figure |4] demonstrates that the force direction as well 
as the magnitude change with t^max ■ Depending on uJmax 
the angle ^p < 9 ov ij) > 9. Since uJmax depends on the 
particle speed v and the distance zq from the dielectric, 
the force magnitude as well as the force direction also 
depend on v and zq. This is in strong contrast to the case 
of a dielectric with no losses and no dispersion, where the 
force direction is independent of zq. 

Interesting behavior is observed when uJmax > i^pi > 
a;p2, fig. [5] The angle ip is negative for low 9, which 
means that the y components of the force and velocity 
have the same sign and the particle is accelerated in the 
y direction. Nevertheless, this does not violate energy 
conservation because the product v • Fp < and the 
particle's energy decreases. 

Peculiar behavior of the transverse force is observed 
for values of 9 below « 20° where it becomes repulsive 
(positive). This is shown in the bottom panel of fig. pi 
where we split the interaction into the Cerenkov part 
(a;p2 < uj < ujpi) and the evanescent part (0 < uj < ujp2 
and UJ > ujpi). For this particular case the Cerenkov con- 
tribution, which can become repulsive, starts to dominate 
for low angles. The origin of the repulsion is the mo- 
mentum carried into the dielectric by the excited waves, 
balanced by the particle [TU]. It is therefore possi- 
ble for a dielectric surface not only to repel a trans- 
verse charge packet as demonstrated in ref . [13] , but also 
a point charge; an outcome which seems to contradict 
common notions based on electrostatic considerations of 
point charges above metallic or dielectric surfaces. 
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